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A scheme to achieve spin squeezing using a geometric phase induced by a single mechanical mode 
is proposed. The analytical and numerical results show that the ultimate degree of spin squeezing 
depends on the parameter , which is the ratio between the thermal excitation, the quality 

factor and square root of ensemble size. The undesired coupling between the spin ensemble and 
the bath can be efficiently suppressed by Bang-Bang control pulses. With high quality factor, the 
ultimate limit of the ideal one-axis twisting spin squeezing can be obtained for an NV ensemble in 
diamond. 
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Introduction - The NV centers in diamond are amongst 
the most promising implementations of quantum bits for 
quantum information processing [l| and nanoscale sen¬ 
sors Q, which is because their ground state spin triplet 
posses ultra-long coherent time at room temperature Q 
and can be readout via optical fluorescence. Significant 
progresses have been achieved in recent experiments to 
couple the NV electronic spins to nuclear spins [ 3 , @ and 
mechanical resonators [3,j3L' The nanoscale magnetom- 
etry ii , thermometer [l^ and electric field detection 


llj have been demonstrated by single NV or an ensem¬ 


ble. 


It’s well known that the quantum states can boost 
the precision measurement beyond the standard quan¬ 
tum limit I12II . Among them, the spin squeezed states 
(SSS) |13l - [l^ have attracted a lot of interest and applied 
to spin or atom ensembles for atomic clocks and gravita¬ 
tional wave interferometers. There are many proposals 
and experiments to realize the spin squeezing in atom 
ensembles, such as atom-atom collisions 11711 . quantum 
non-demolitionfQND) measurement 0 llflj l and cavity 
squeezing [20l - [^ . Very recently, spin squeezing of an NV 
ensemble by Tavis-Gummings type interaction between 
phonon and spins has been proposed for quantum 
enhanced magnetometry. 

In this paper, we propose a new approach for the re¬ 
alization of spin squeezing by phonon induced geomet¬ 
ric phase, using an ensemble of NV centers dispersively 
coupled to a mechanical resonator. It’s shown that the 
ultimate degree of spin squeezing by one-axis twisting 
can be realized, for reasonable ratio between the thermal 
excitation and the quality factor of mechanical oscilla¬ 
tors. Furthermore, the effect of the coupling between 
NV centers and environment is studied, which leads to 
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dephasing and degrades the spin squeezing effect. By in¬ 
troducing Bang-Bang pulses sequence, the decoherence is 
effectively suppressed and significant spin squeezing can 
be achieved for the NV ensemble. 

Model- The negatively charged NV center (NV“) 
in diamond is well-studied, whose Hamiltonian reads 

Hnv = {D+S\e^)Sl+tibgeS P,!!!, where D « 2.87 
GHz is zero-field splitting, and Cz are axial ground- 
state electric dipole moment and electric field (strain 
field), respectively. With appropriate bias field Bz, the 
two microwave transitions |0) o |±1) can be addressed 
separately in experiment, and we focus on the |0) O |—1), 
which transition can be treated as a spin-^ system in the 
following. Putting the NV“ spin ensemble in a gradi¬ 
ent magnetic field ^ 0, then the displacement of 
diamond or nanomagnet 5u will shift the transition fre¬ 
quency by Aujnv = ThgeCTz^^du [2^. Alternatively, the 
strain field of a diamond nanomechanical oscillator can 
induce an electric field inside the cr yst al and giye rise to a 
similar phonon-spin interaction |29| . Both approaches 
to couple the spin with nanomechanical oscillator haye 
been demonstrated in experiments recently [3,0, . 

The simplified Hamiltonian of an ensemble of 2V spins 
coupled to a mechanical resonator is 

H = uJaa'''a-\-gJz {a-\-a^) , (1) 

where a and are annihilation and creation operators of 
phonon, uia is the frequency of the mechanical resonator, 
Jz = ^ collectiye spin operator, and g 

is the single phonon coupling strength. Along with the 
progresses in the nanofabrication of diamond material, 
yarious diamond nanomechanical resonators haye been 
realized in experiment, with frequency ranging from 1 
kHz to 1 GHz, and the quality factor Q ranging from 
100 to around 10® [1, [13, l33 - l33| . Thus, we study the 
spin squeezing induced by the mechanical resonator with 
frequency ujg /27r = 1 MHz and coupling strength g/2TT = 
1 kHz 13511 in this work. 
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FIG. 1: Trajectories on phase space of a coherent wave packet 
(a)=Re(a) + ilm(a) for spin state |m) with m — 0, ±1, ± 
2, ± 3. Here we set g/oJa as a unit. 


The Hamiltonian preserves Jz of the spin ensemble. 
When we integrate the Schrodinger equation for collec¬ 
tive spin states, the phonon will introduces a spin de¬ 
pendent geometric phase shift [s^. It is convenient to 
study the mechanical resonator by the coherent state 
|a), and we can write ja) = |Re(Q;)-I-Hm(Qf)) with 
a{t) = [l —The coherent state behaves 

somehow iflce classical particles in phase space. Its cen¬ 
ter, given by Re(a) and Im(a), follows a classical tra¬ 
jectory, while the width of these wave packets remains 
fixed, which is given by the uncertainty of the Re(a) and 
Im(a). In Fig. I we plot the usual phase-space trajec¬ 
tories for (a) = Re(a) -I- Hm(a), and we have used the 
eigenstates \m) of spin operator Jz as the initial states 
of the collective spin. We plot phase-space trajectories 
only with m = 0, ±1, ±2, ±3 for simple explanation, 
which clearly show that the coherent wave packet is re¬ 
stored to its original state after a fixed time T = ‘I'k a 
or integer times of T. For different |m|, there are differ¬ 
ent radius circles in the phase-space trajectories, and it 
is the central symmetry for the opposite m. 

The geometric phase, as the enclosed circle area of the 
traj ectory in phase space, is insensitive to the initial state 
[^. Thus, phonon induced geometry phase is robust 
against the imperfection of initial phonon state prepara¬ 
tion, and we assume a (0) = 0 for simplicity. However, 
the decay and thermal noise of phonon during the spin- 
phonon interaction will influence the geometry phase ac¬ 
cumulation. In this case, the system dynamics follows 
the Master equation 

^ [H, P] + I {nth + 1) C{a)p + ^nthC{a))p- (2) 

Here ^ = lo a!Q describes the decay rate of the mechan¬ 
ical mode, rith is the mean phonon number of the me¬ 
chanical thermal noise and C{o)p = 2opo' — 6^op — po^o 
is the Lindblad superoperator for given jump operator o. 
The reduced density matrix of the collective spin can be 
written as 

Psprn = |to) {n\ . (3) 

m,n 


4'm,n{t) is the phase difference between these spin states. 
The phase can be solved as 


4‘m,n{i) 

= - [nth + 0 “ nif | a(T) dT+ \ a{t) |^| 


+ig{n^ 


— TO^)Re / a{T)dT. 

Jo 


( 4 ) 


Here, the amplitude of mechanical resonator is a{t) = 
■y/ 2 +iui [4 ~ 371. The finite 7 of the me¬ 

chanical resonator introduces decoherence and leads to 
the first term of the above equation, the second term is 
corresponding to the interaction inducing spin squeez¬ 
ing. Molmer and Sorensen proposed an approach for 
ion-trap to realize the spin squeezing, which is insensi¬ 
tive to the initial thermal phonon states [1^. Compare 
to the Molmer-Sorensen scheme that two laser pumping 
and Lamb-Dicke approximation are required , our ap¬ 
proach utilize stable spin-phonon interaction and there is 
no approximation in our model. 

Spin squeezing.- The spin squeezing is evaluated by 
squeezing parameter [H, [l 3 


^2 niin{^Jl^) 


( 5 ) 


where A is the variance of spin operators along direc¬ 
tion perpendicular to the mean-spin direction rio = , 

which is determined by the expectation values (Ja), with 
a G {a;,y,z}. For an atomic system initialized in a co¬ 
herent spin state (CSS) [S^ along the x axis, satisfy¬ 
ing Jx l ^(0))atom = IV-(0))atom ^ ^e have Pm,n (0) = 

q-2Ar / (2N)1 (2N)! , A 7^ — /V/9 

^ yw — m)\{N-\-m)\ {N—n)\{N-\-n)\ fi± / * 

Thus, for squeezed spin states we have < 1. 

First of all, we studied the spin squeezing by Eq. (3) 
without thermal noise. The squeezing parameters as a 
function of the time (dimensional number gt) for various 
quality factor Q are plotted in Fig. 2(a). As expected, 
the effect of phonon induced geometry phase leads to the 
twisting and squeezing of CSS, thus the decreasing 
with time. After a certain optimal t, the increases, 
due to the over twisting by the geometry phase, and high 
order effect arises. It is shown that the minimal value of 
spin squeezing parameter decreases with higher mechan¬ 
ical quality factor Q. When the quality factor Q = 1000 
(the black solid line), the almost perfect spin squeezing 
for the ideal one-axis twisting can be achieved. Including 
the mechanical thermal noise, squeezing parameters 
as functions of the time with the quality factor Q — 1000 
are plotted in Fig. 2(b). It is natural that the spin 
squeezing becomes worse with the increasing of the ther¬ 
mal noise nth- We also studied the suppression of the 
influence of thermal noise by improving the quality fac¬ 
tor Q. As shown in Fig. 2(c), the optimal spin squeezing 
(the minimum value of the ^%{t)) is plotted against the 
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FIG. 2: (Color online) (a) The squeezing parameter as a 
function of the time for nth = 0 and various Q = 5 , 10 , 
1000 (top to bottom), (b) The squeezing parameter as a 
function of the time for Q = 1000 and various nth = 100 , 
50 , 0 (top to bottom), (c) The green solid line is optimal 
squeezing parameter versus the quality factor Q for nth = 
100, and the red dashed line is the result for ideal one-axis 
twisting spin squeezing. = 10 for all simulations. 


Q for nth = 100. The reduces with Q and approaches 
the limit for ideal one-axis twisting spin squeezing (red 
dashed line). 

To understand these results, we simplified the spin 
state state dependent geometric phase 


0 m,n(t) — 


(7/2)^ -bOJ 


(m^ — n^) — nf , ( 6 ) 


under the approximation t 7 which means a{t) = 
— and the transient evolution of the mechanical 
resonator is neglected. Here, the dimensionless factor 
/i = iuh+ 1 / 2 ^ accounts for the coefficient 

proportional to the time t, and the two terms within 
the bracket corresponds to spin squeezing and decoher¬ 
ence, respectively. Then, we can obtain the degree of 
spin squeezing for the initial state CSS along the x axis 

= (7) 


squeezing scales with From Eq. (7), we obtain the 

approximated upper bound of the optimal spin squeezing 
^ ^/(l —e ^), which indicating that 

the ratio should be as small as possible. 

As long as < 10“^, we can achieve squeezing al¬ 

most as good as the best squeezing achievable with ideal 
single axes twisting [Fig. 2(c)]. 

Bang-Bang control. - During the preparation of optimal 
SSS, there are inevitable couplings between the system 
and baths. For example, the lattice vibrations and envi¬ 
ronment spins will induce dephasing and destroy the spin 
squeezing. The dynamical decoupling technique is well 
known for protecting coherence from environment koj- 
Id^ . and now we apply the Bang-Bang (BB) pulses |4Cll | 
to suppress the decoherence. The sequence consists of M 
pulses, which split the total time interval t into M small 
intervals tp = with p = 1,2,...,M. The pulses rotate 
the collective spin states around y axis, and we chose 
the pulse sequence to rotate tt and —tt alternately, which 
leads to = —azj and = j 2 ^ 

Therefore, the spin squeezing is conserved while the az 
is inverted by the BB pulses. Considering the 2N qubits 
which are independently coupled to thermal baths, the 
Hamiltonian from the Eq. o is changed to 

H' = ojaa^g -\- ge{t)Jz {a + Q^) + i^kb\bk 

k 

2N . . 

+ — ^^fT-kjih + bl). (9) 

j=l k 

Here, bk and b\. are the creation and annihilation bosonic 
operators of the fc-th bath mode, which coupling to the j- 
th spin with coupling strength h^j. The switch function 
e(r) due to BB pulses is given by e(r) = 
tp)d{tp+i — t) with 9 {t) is the Heaviside step function. 

With the decoherence and BB, the geometry phase fac¬ 
tor [Eq. ([3])] of spin states are solved as 

= - (^nth + 0 I 1^ dT+ I a'{t) p| 

-\-ig{n‘^ — m‘^)Re f e{T)a'(T)dT — Km,n(t). (10) 

Jo 


where 

A = l-cos^^-'^{2Ct) 

B = -4sin(C't)cos2^-2 (Ct) ( 8 ) 


Here, a'{t) = ~ig and K.m,n{t) 

is due to the decoherence. Assume that the baths to 
each spin are Ohmic and have the same spectral density 
we have 


Here, Ct = ^ 17774 ^ x gt- The analytical solution im¬ 
plies that the spin squeezing is mainly determined by the 
two dimensionless parameters Ct and p,. For Q ^ 1., we 
have 4 _)_ 4 y 4 Qa ~ 1. For N ^ 1, we can apply the approx¬ 
imation cos^^“^ (x) Ki for a: <C 1. So, the time 

required {gt ~ 160 in Fig. 2) to achieve the optimal spin 


Km,n(t) < (|n - to]- b 2 ) / G (oj) FM{oJ,t)duJ, (11) 
Jo 

where the modulation spectrum is FM{w,t) = 
( 2 M+ 2 )(i+^( 1 ) — cos(qjt)) ^ temperature-dependent 

interacting spectrum is G(w) = riu)e~~ ( 1) 
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FIG. 3: (Color online) (a) The squeezing parameter as a 
function of the time gt. The parameters are g — 0, M = 0 
(black line), ?7 = 4x 10“"^, M — 0 (red line), and g = 4x 10“'*, 
M = 500 (blue line), (b) The blue solid line for the optimal 
squeezing parameter versus the pulses M for = 4 x 10“*, 
and the black dashed line for the result without thermal baths. 
Other parameters are nth = 10, \ = 4/g, Q = 1000 and 
N = 10. 

, in which g is the coupling strength between the sys¬ 
tem and the bath modes, lUc is the cutoff frequency, and 
A = 1/kbT is the inverse temperature. In order to sim¬ 
plify the calculation, we use the upper limit instead of 

the Km,n(t). 

In Fig. 3(a), we numerically calculated squeezing pa¬ 
rameter as a function of time for various g and M. 
Since the decoherence term Km,nit) is proportional to 
the coupling strength between the system and the bath 
modes, we observe the incremental of the optimal 
for increasing g (black and red lines). The blue line 
shows the suppression of decoherence by BB, and here we 
choose the sequence number M = 500 and g = Ax 10“*, 
which is contrast to the red line. There are periodic peaks 
with the separation distance Agt = tt, and the peak val¬ 
ues are obtained when gt = (n-|-l/2)7r, n is integer. 
This phenomena can be interpreted as following: the BB 


pulse period is Tm = t/M, and the time period for the 
phonon state trajectories in the phase-space [Fig. 1] is 
T = For M = 500 and uja/g = 1000, we have 

Tm/T — gt/TT. When gt/i: = n is integer, the geometric 
phase is always cumulative, and the coherent spin squeez¬ 
ing effect is not degraded by the BB pulse sequence. In 
contrast, when gt/ir = n + ^, the geometric phase im¬ 
prints alternating sign as function of M , and then the 
spin squeezing is weakened. Compare the minimas of 
with BB (blue line) to the results without BB (black 
and red lines), the undesired effect of decoherence is effec¬ 
tively suppressed by the dynamical decoupling. Fig. 3(c) 
shows the optimal versus the pulse sequence length 
M. With increasing M, the optimal value of is im¬ 
proved and approaches the red dashed line, which is the 
ideal result determined by the Eq ([51) without thermal 
baths. When M > 400, the influence of thermal baths 
on spin squeezing can be almost eliminated, which means 
Kni,n (t) ~ 0. 

Conclusion.- An approach to achieve spin squeezing 
by phonon induced geometric phase is proposed. This 
scheme is feasible for experiments on solid state spin en¬ 
semble coupled to a mechanical oscillator. With reason¬ 
able parameters, the ultimate limit of the ideal one-axis 
twisting spin squeezing can be achieved as long as the 
quality factor is sufficiently high that Q > x 10^. 

The decoherence due to spin-bath coupling can be ef¬ 
fectively suppressed by the Bang-Bang pulses. This 
geometric-phase-based spin squeezing can be used to sig¬ 
nificantly improve the sensitivity of magnetic sensing 
with nitrogen-vacancy spin ensembles. Moreover, the 
technique can be generalized to spin ensembles coupled 
to other high-Q Bosonic modes that prepare quantum 
states by geometry phase. 
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